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This paper reviews the status of Lattice QCD where relevant to kaon physics. Some summary 
intended for non-experts of the broader status of the field is given in section 1 . I interpret this goal 
as including some aspects of pion physics and the Chiral Lagrangian in section 2. 

I review the status of Lattice QCD input to determining Vus, section 3, both via fx/ fn, sec- 
tion 3.1, and Kij, section 3.2. The kaon bag parameter Bk which contributes to the Ek constraint 
on the unitarity triangle is reviewed, section 4. Lattice K — nn calculations were reviewed at this 
conference by Norman Christ [1] and are thus beyond the scope of this review. 

Carrying out a review of this nature is an increasingly daunting task in the presence of the 
excellent Flavia Lattice Averaging Group effort [2] . It is also worth mentioning excellent recent 
reviews by individual members of FLAG [3, 4]. In many areas it is difficult to match this com- 
prehensive effort. My approach therefore is to try to highlight theoretical issues, distinguishing 
features, and methods of the lattice calculations that are not easily conveyed by the compressed 
asterisk rating system; I do not necessarily attempt the comprehensive categorisation and averag- 
ing task this working group has undertaken except where I may have significant updates from new 
material, or wish to differ in the analysis. 

PAB wishes to thank the many authors whose original research has been reviewed and empha- 
sizes that any use of Lattice world averages or quoted results herein must cite the orginal source(s) 
where the intellectual value was created. 

1. Current Large Scale Lattice Simulations 

Despite many recent advances in both algorithms and computing hardware Lattice QCD can- 
not presently be performed at physical masses in a large volume for reasons of numerical cost. 
Were this feasible, there would be little debate about methods in the field. In practice, the choice of 
approach, table 1 , is a complex optimisation problem where the selected approach is dependent on 
the physical processes for which one aims to yield a minimal overall eiTor. These eiTors are com- 
posed of statistical and finite volume, mass extrapolation and discretisation systematics. Strategy is 
influenced by estimates of the importance & calculability of finite volume, NNLO chiral expansion 
effects, and even the influence of strange quark loops. 

The physics goal of computing chirally structured weak matrix elements (Bk and K Tin, for 
example) makes perhaps the biggest rationally justifiable difference of approach. RBC-UKQCD 
and JLQCD use substantially more expensive methods to obtain near exact chiral symmetry and 
gain access to richer phenomenology than would otherwise be possible. These approaches also give 
automatic 0{a) improvement off-shell, making improved Rl-mom non-perturbative renormalisa- 
tion simpler. However, less theoretically pristine approaches are substantially cheaper and BMW 
and MILC in particular have used this to gain access to finer lattice spacings and lighter simulated 
masses in large volumes. 

Lattice QCD gives complete freedom in choice of the quark masses (dynamical loops and va- 
lence legs) in numerical simulation. Purely unitary simulations maintain the valence quark masses 
equal to those of simulated dynamical flavors in all matrix elements. A partially quenched simu- 
lation is a superset of the corresponding unitary simulation; here the multiple valence masses ai^e 
used (both equal to and not equal to the dynamical masses) and next-to-leading order PQ-;^''^is 
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Coll3bor3,tion 


Action 


L(fm) 


m™'" (MeV) 


fl-i GeV 


ChPT 


Unitary 


PV corr. 


JLQCD[5] 


Overlap 


1.7 


310 


1.8 


SU(2) 


U 




CDH 


RBC-UKQCD [6] 


DWF 


2.0,2.7 


330 


1.7,2.3 


SU(2),SU(3) 


PQ 






ETMC [7] 


Twisted mass 


2.0 


300 


1.9,2.2,2.8 


SU(2) 


U 




CDH 


PACS-CS[8] 


Clover 


3.1 


156 


2.2 


SU(2),SU(3) 


U 




CDH 


BMW [15] 


Clover 


4.0 


190 


1.6,2.3,3.03 


polynomial 








MILC[16] 


Staggered 


2.4,2.9 
3.4 


320 


1.1,1.3,1.6 
2.2,3.3,4.4 


rs-SU(3) 


PQ 







Table 1: Summary of the parameters of major lattice calculations reviewed 



available for a variety of important quantities such as pseudoscalar masses, decay constants and the 
kaon bag parameter. 

RBC-UKQCD and MILC have used partially quenched analyses to increase the information 
within the chiral regime, while JLQCD, ETMC, BMW and PACS-CS used unitary simulations. 

Where smaller volumes are used, it is common for finite volume corrections to be applied to 
the data prior to chiral extrapolation. This, of course, takes the effect of finite volumes as an input 
to (rather than output of) the lattice calculation. The presently favoured model is the resummed 
scheme of Colangelo, Durr and Haefeli (CDH) and is applicable only to unitary datapoints, and 
cannot be combined with partially quenched analyses. For this reason the calculations where mj^L 
becomes rather smaller than around 3.2 both apply finite volume correction and consider only 
unitary datapoints. 

Approaches to mass extrapolations include the application of SU(3) and SU(2) chiral effective 
theories. BMW recently adopted Taylor expansion having achieved both light simulation masses 
and a large volume. They feel sufficiently close to the physical point to include polynomial fits in 
their analysis; these ai^e viewed as an analytic expansion around the non-zero physical quark mass 
in a massive region whose validity does not extend to the chiral non-analyticity at zero. 

1.1 Is lattice QCD ab initio in practice? 

Lattice QCD has rightly enjoyed the reputation of an ab initio method to solve QCD. However, 
present analysis & extrapolation approaches do involve varying degrees of «o«-lattice input. What 
is clear is that, in some cases, the extrapolated results have relinquished the significant benefit of 
having only the QCD Lagrangian as input. 

Ideally, any use of the low energy expansion of QCD for chiral extrapolation should be per- 
formed in a consistent & well motivated context. Either complete NLO or complete NNLO should 
be used; it should both describe the data and appear a well convergent series within the range of 
the data without ad hoc resummation schemes or partial NNLO where only favoured analytic terms 
are included, or certain EEC's biased towards phenomenological input. 

Taking the long view, lattice QCD remains systematically improvable. Continued advance in 
supercomputing power will soon remove the need for such extrapolations, just as the quenched 
approximation has been long discarded. However, it is worth highlighting mass extrapolation and 
finite volume issues as areas where the field should do better. 
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2. Chiral Effective Lagrangian 

The chiral effective theory can be viewed in two roles by a lattice practitioner. First it can 
be viewed as an important predictive tool which can be verified and understood with lattice QCD 
(i.e. intrinsic interest). The chiral effective theory has also been used as a critical component of 
making well founded mass extrapolations from simulation in almost all recent lattice calculations 
(i.e. practical value). The former requires that chiral perturbation theory be a well convergent series 
at the relevant physical masses, while the latter requires convergence above the physical masses. 

Were SU(3) chiral effective theory convergent at the kaon mass, it would be the natural frame- 
work for interpreting lattice data. MILC in particular have relied on the partially quenched SU(3) 
approach to remove unitarity violating discretisation effects introduced by their "rooting" prescrip- 
tion [17]. However, some determinations of the leading order low energy constant Fq representing 
the pseudo-scalar decay constant in the SU(3) chiral limit (table 2) have been surprisingly low. 

RBC-UKQCD found [6] that NLO SU(3) formula does not describe pseudoscalar masses and 
decay constants above niji ~ 400MeV. In addition to surprisingly low results for the LEC, RBC- 
UKQCD obtained a slightly low result = 124.1(3.6)(6.9) MeV ; with a single lattice spacing 
this is within systematic error of the physical value. At simulated masses and even at the pion 
mass the SU(3) correction to leading order was surprisingly large, figure 1. For this reason, they 
introduced the appUcation of SU (2) ;f ^^to both pionic and kaonic quantities in an underlying 2+ 1 
flavor simulation [11,6], and this has been adopted by other groups. 

In this SU(2) effective theory the kaon is treated as a heavy external meson field [11, 12, 6]. 

The least convergent terms expand in ( ) , and this appears the more convergent approach for 



The R\CS-CS collaboration also found that kaon data was better described by SU(2)[13, 14], 
figure 2; however their result for Fq was much larger and the conclusion about simulated 200- 
400 MeV pions differs. The underlying data for is qualitatively similar between R\CS-CS and 
RBC-UKQCD but differs by around 4%. R\CS-CS determine Za from one loop lattice perturbation 
theory; whether both the overall scale and curvature of their chiral fits would become closer with a 
non-perturbative Za remains an open question. 

The MILC collaboration found a larger Fq than RBC-UKQCD [16]; however this is not nec- 
essarily discrepant when one considers that the RBC-UKQCD results were from a single lattice 
spacing. The SU(3) chiral limit could be clarified relatively easily by simulations with artificially 
light "strange" quarks.. 

2.1 The curious linearity of lattice data 

BMW have not yet published a detailed comparison of their analytic fits to their chiral fits 
[15]. However several collaborations have found suprising linearity in their data for f^. 

Since the Kaon conference RBC-UKQCD have presented a second lattice spacing [18, 19] 
at Lattice 2009, where they found good scaling compared with their earlier calculation and form 
conclusions about the continuum limit with theoretically pure dynamical chiral fermions. They set 
the scale using the Q. baryon. Figure 1 demonstrates that while a simple linear- from their 290- 
450MeV simulated masses produces the PDG fji, NLO extrapolation does not. The discrepancy is 




current data than SU (3) expansion where the least convergent terms expand in 
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Collaboration 


Fo (MeV) 


F/Fo 


MILC [16] 


106(8) 


1.15(5)(+J3^) 


RBC/UKQCD [6] 


93.5(7.3) 


1.229(59) 


PACS-CS[8] 


118.5(9) 


1.078(58) 



Table 2: There remains some spread in recent results for the leading order low energy constant Fq of the 
SU(3) chiral effective theory. A low value would imply SU(3) corrections for kaons in particular must be 
uncomfortably large. F is the SU(2) low energy constant and here the normalisation is such that the pion 
decay constant is around 131 MeV. 

consistent with a naturally sized NNLO effect; however the cancellation to produce linearity within 
statistical error is striking. This data corroborates the BMW claim that Taylor expansion in the 
region directly above the physical masses shows good convergence; however the presence of known 
non-analytic terms makes it difficult to continue analytic fits to lighter masses and RBC-UKQCD 
used the difference between analytic and chiral fits as a systematic error in predictions [18]. 




140 250 330 420 nj^p'-" 



Figure 1: Left: NLO chiral extrapolation of RBC-UKQCD's dynamical DWF decay constant for fj^ 
on their 1 .73GeV ensemble. The decay constant in the SU(3) limit is surprisingly low, and corrections in 
this expansion are uncomfortably large. This plot displays partially quenched data. Right: Since the Kaon 
conference RBC-UKQCD have presented a preliminary continuum limit for their simulations. The data 
points here are adjusted to represent values in the continuum limit and compared to both an NLO and linear 
extrapolation. The scale is taken from Mq, and NLO j^^yields a pion decay constant that is inconsistent with 
its physical value, while linear extrapolation is in agreement. The inconsistency is certainly explicable by a 
naturally sized NNLO effect, however the coincidence of cancellation between orders to leave such linear 
behaviour at this level of precision is curious. This plot displays only the unitary datapoints for simplicity. 

2.2 Low energy constants 

Lattice determinations of the leading order and next-to-leading order LEC's have been sur- 
veyed in detail by S. Necco, and by FLAG in excellent recent reviews [4, 2]. I reproduce here 
some of the SU(2) NLO summary in table 3, as I believe this is more convergent and hence more 
meaningful than SU(3). I recommend arXiv:0901.4257 [4] for details and further results, including 
a survey of the large body of e-regime results which ai^e not covered in this paper. 
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Figure 2: PACS-CS NLO fits to the kaon decay constant using SU(2) (left) and SU(3) (right) . As with 
RBC-UKQCD SU(2) is found to give better agreement at larger masses. The vertical dashed line represents 
the physical point, albeit in a small volume. The "FSE" fits are corrected for finite size effects using CDH. 
SCALE? 
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Collaboration 
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h (SU(2)) 


U (SU(2)) 


'"3 (SU(3)) 


U (SU(3)) 


ETM 


2 


3.42(8)(10)(27) 


4.59(4)(2)(13) 






ETM 


2 


3.2(4)(2) 


4.4(1)(1) 






JLQCD/TWQCD 


2 


3-44(57)(+«)(+Jo2) 


4. 14(26) 


+32-J 




MILC 


2+1 






i.i(6)(+;:« 


4.4(4) 


RBC/UKQCD 


2+1 


3.13(33)(24) 


4.43(14)(77) 


2.87(28) 


4.10(5) 


PACS-CS 


2+1 


3.14(23) 


4.04(19) 


3.47(11) 


4.21(11) 



Table 3: (Abridged) summary tables from [4], which also reviews SU(3) LEC's and the large body of 
work making use of the e-regime. Here F is in the convention where fji ~ 92MeV 



2.3 Pion form factors and h 

The pion vector form factor is related to the low energy constant 

(n*{p')\V,\n*{p)) = F"{q^)(p,+p'^) 

A summary of recent results for the pion charge radius is given in table 4. There are presently 
considerable discrepancies in the results at intermediate masses prior to chiral extrapolation, and 
consequently differing conclusions. 
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ETMC and JLQCD [24-27] obtain somewhat lower values than RBC-UKQCD and QCDSF 
in the region 300 < rtij^ < 500 MeV, and conclude the NNLO effects in the chiral expansion ai^e 
significant. This calculation typically involves the use of a number discrete Fourier modes for q 
with model dependent interpolation used to obtain the derivative at = 0. RBC-UKQCD and 
ETMC [30, 24] obtained good resolution in the low q^ region, figure 3. Here, a twisted boundary 
condition with phase e^^ applied to valence Fermions on the simulated torus allows to vary lat- 
tice momenta smoothly between the standard periodic and anti-periodic Fourier modes, enabling 
several momenta between the zero and first Fourier modes [28]. RBC-UKQCD found that NLO 
;^^^yielded good results for the pion charge radius, but this remains a puzzle however, since ETMC 
also used twisted boundary conditions but concluded NNLO ;^^^was required. 

The determination of the pion scalar charge radius by JLQCD is worth noting and which has 
been performed including the effects of disconnected quark flow diagrams for the first time[26]. 



CoUab 


Action 




{rl)s 


ETMC [24] 


2f TM 


0.456(38) 




JLQCD [25, 26] 


2f Overlap 


0.409(23)(37) 


0.617(79)(66) 


JLQCD [27] 


2f Clover 


0.396(10) 


0.60(15) 


QCDSF [29] 


2f Clover 


0.441(19) 




RBC-UKQCD[30] 


2+lf DWF 


0.418(31) 




LHPC [31] 


2+lf MILC/DWF 


0.310(46) 





Table 4: Summary of recent results for the pion charge radius from lattice QCD 
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Figure 3: LEFT: BMW's chiral extrapolation of fK/ fn [15]. This includes an impressive 190MeV lightest 
simulated mass, 4 fm volumes and a continuum limit. Detailed exposition in print is absent, but nevertheless 
this is an exciting calculation for which further details are eagerly awaited. RIGHT: RBC-UKQCD mapped 
out the deep infra-red region of the pion form factor using twisted boundary conditions. ETMC used a 
similar approach down O.OSGeV^ [30] 



2.4 L\, 



Here we see the first of many quantities that ai^e made accessible with chirally symmetric lattice 
actions. L'jg is related to the electromagnetic component of the splitting between and n^, and is 
determined from small difference between the axial and vector vacuum polarisation functions. 
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V-A 



^V-A--f-^L''M 24712 

JLQCD produced a beautiful first determination in their «y = 2 overlap simulation [32]. obtaining 
the pion e-m mass squared splitting (in chiral limit; c.f. PDG with physical u^d masses 1261MeV2): 

^^oK) = -5-2(2)1^ X 10-3 

»4± - m^o = 993(12)(;?35)(149)MeV2 

A second calculation has been performed by RBC-UKQCD with 2+lf domain wall fermions 
[33]. The calculation is also interesting in that it produces information about the radius of conver- 
gence of x^in momentum space, in addition to its convergence in pion mass; these are similar. 



The two best lattice constraints for V,,, involve SU(3) breaking effects: for ^ this is 0(20%), 
and for f§^{q^ = 0) this is 0(4%). While the Kit, approach is less mature, it looks very promising 
and is rapidly becoming better studied. We consider these two quantities in turn. 

3.1 Vus^romfx/fn 

This topic was reviewed in detail by Lellouch [3]. The calculation with the smallest quoted 
eiTors is the mixed action HISQ/Asqtad staggered calculation by HPQCD[35], followed by the 
smeared clover simulation from BMW[15]. 

The HPQCD calculation used a mixed action. While published in an abbreviated form, it ap- 
pears to have not followed MILC's best practice of performing a partially quenched rs-;^^''^ analysis 
to absorb unitarity violations. Rather, it appears to have matched the HISQ valence pseudoscalar 
Goldstone taste in to the (irrelevant here) AsqTad valence pseudoscalar Goldstone taste. This dif- 
fers from the various relevant masses of sea pseudoscalar multiplet. It is hard for lattice theorists 
not directly involved in the calculation for form a clear judgement with the available information. 

Figure 3 (left) reproduces the truly impressive BMW simulation that spans three lattice spac- 
ings and includes volumes up to 4 fm and pion masses down to an impressive 190MeV. The 
BMW calculation has not yet received a detailed exposition in print, and has so far lacked a non- 
perturbative determination of the axial cun^ent renormalisation so important checks on fK and fn 
seperately have not yet been performed; it seems likely non-perturbative operator improvement 
was not used and resulted in their dual use of a and continuum extrapolation. Their chiral ex- 
trapolation was interesting but details ai^e sparse; chiral expansion and Taylor expansion fits were 
combined weighted by quality of fit and only the combined fit has been shown. Publication of the 
qualities of fit, and direct comparison of the two in the style of figure 1, would be interesting since 
continuing an analytic fit downwards is not without risk and it is not yet clear how much the lightest 
data points constrain chiral curvature. 
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Since the Lellouch review, ETMC have updated [36] their Nf = 2 twisted mass resuU to 
Ik/ fn = 1.210(18). As this will not substantially change the lattice average I recommend con- 
tinued use of that quoted by Lellouch [3] while noting the leading publications are in formats that 
are not amenable to critical assessment 



A/A = 1.194(3) (10). 

3.2 

The semi-leptonic form factor 

W)mK{p))=U{q^){p^+p'^)+nq'){p^-p'^) 
can be computed precisely using several "double ratios" following [37, 38], such as 
{K{Q)\Vo 1 7r(0) ) (ii:(0) I Vol 71(0) ) [niK + m^Y 



(^(0) I Vo 1^(0) ) (tt (0) I Vo 1 71 (0) ) "^mKin^ 



-\h{qL.)Y 



Ref. 



Nf action a;[fm] L[fm] 



typ/val 



/+(0) 



JLQCD[40] 2 NPSW 0.09 1.8 

RBC[41] 2 DWF 0.12 2.5 

ETMC [42] 2 tmQCD 0.11 2.7 



FNAL/MILC[44] 2+1 KS+Wil 
RBC/UKQCD[43] 2+1 DWF 



0.11 1.8, 2.e 



550/550 0.967(6) 
490/490 0.968(9)(6) 
260/260 0.958 1(57)(35) 

0.962(6)(9) 
290/240 0.9644(33)(34)(14) 



Table 5: An update of the corresponding table by Lellouch [3] to include new ETMC results 




Figure 4: LEFT: RBC-UKQCD chiral extrapolation of /+(0) represents the best lattice constraint and has 
2+1 flavors. RIGHT: Courtesy A. Jiittner - using RBC-UKQCD's result for Fq to better define /2 results in 
NLO describing the data well and the chiral expansion of this small breaking looking naturally convergent. 
A redefined would be near zero. 
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The form factor f+{q^ = 0) is of key interest for Vus and has chiral expansion 1 +/2 +/4 • • • 
where /„ ~ 0{M'^ j^/ {AtzFqY). The best constraint in, table 5 uses 2+1 flavors and is from RBC- 
UKQCD [43]. Their chiral extrapolation is shown in figure 4 obtaining 

/+(0) = 0.9644(33)''^'(34)'='"^^P°''^^^™(14)''''^ 

The calculation used Fourier modes, with model dependence in interpolation and chiral extrapo- 
lation; there was no controlled continuum limit and a discretisation systematic error was budgeted 
at 4% of 1 — /+ (0) . The historical use of fji as the denominator for /2 is not the only sensible 
choice. Taking the above RBC-UKQCD result for Fq in table 2 better matches the data and re- 
arranges the expansion at NNLO and above. Different reasonable choices of /2 seem to be an 
important systematic in SU(3) based extrapolations of the form factor. RBC-UKQCD presented 
additional data at the lattice conference making use of twisted boundary conditions to simulate 
directly at = [39, 45] and have developed an SU(2) formalism [46] which they will use in the 
chiral extrapolation of their new data. 

The ETMC fif = 2 calculation [42] had a notably robust chiral extrapolation with pion masses 
as low as 260 MeV and made use of both SU(3) and SU(2) formalisms in their chiral extrapo- 
lation. An adjustment was made (but not quoted here) to "coiTcct" the lattice calculation for the 
missing strange quark using the SU(3) expression for The degree to which leading effects of 
the absent strange loops are akeady reabsorbed when the lattice spacing is determined makes this 
unconvincing; I prefer to quote the unadjusted ETMC calculation as an excellent fif = 2 result. 

4. Bk 

The matrix element 

l{K^\Ao){Ao\K^) 

is multiplicatively renormalised in the continuum and for lattice actions with chiral symmetry. Non- 
symmetric actions must deal with unphysical taste/chirality mixings which inflate errors. In some 
cases mixed actions are used (e.g. ALV, ETMC), where a valence quark action with multiplica- 
tive operator renormalisation is combined with a cheaper action for sea quarks; resultant unitarity 
violations should be fitted away, typically using mixed action chiral perturbation theory. 

4.1 Bk in the quenched approximation 

For many years the benchmark quenched Bk calculation was from JLQCD using many lattice 
spacings with the staggered fermion action [47]. This calculation involved the lattice perturbative 
treatment of a taste mixed operator basis. The continuum limit of 

B^{nf = 0,2GeV) = 0.565(4) (5), 

obtained by JLQCD [48] using the chirally symmetric domain wall fermion formulation, gives con- 
tinuum limit consistent with various earlier calculations. This includes RBC who used a different 
gauge action thus displaying universality. JLQCD's results are particularly clean in that they apply 
non-perturbative step scaling to enable matching to continuum perturbation theory at large scales. 
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This is substantially below the previous quenched "benchmark", and resolved a puzzle high- 
lighted by Lellouch. It also highlights that lattice renormalisation issues of various formulations 
can be a very important effect. 
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Figure 5: LEFT: Quenched continuum limit of Bk has been substantially clarified, with the previous bench- 
mark using staggered fermions appearing suspect. The chirally symmetric DWF formulation has only mul- 
tiplicative renormalisation and CP-PACS have taken a high quality continuum limit with non-perturbative 
renormalisation. This is plausibly consistent with both RBC (non-perturbative) and earlier (perturb atively 
renormalised) CP-PACS results. RIGHT: RBC-UKQCD chiral extrapolation of Bk in the continuum limit; 
data points are shifted to the continuum limit and overlayed. Systematic error taken from difference between 
linear and SU(2) NLO fits. 



4.2 HBk 

There are two recent and high quality determinations of Bk in the two flavour theory. ILQCD 
[49] made an impressive calculation with 2f overlap action, while ETMC used a somewhat less 
aesthetic mixed action approach with Osterwalder-Seiler valence fermions with 2f twisted mass 
sea fermions [50]. JLQCD obtained B^{nf = 2,2GeV) = 0.537(4) (40) and ETMC obtained 
B^^{nf = 2,2GeV) = 0.56(2) both with single lattice spacings. The ETMC calculaton is prelimi- 
nary and has not yet addressed important unitarity violation systematics and I quote: 

B^{nf = 2,2GeV) = 0.537(4)(40). 

4.3 I+HBk 

Two 2-i-If calculations were considered in detail at the Kaon conference: RBC-UKQCD's 
dynamical domain wall simulation, and the Aubin/Laiho/Van de Water mixed DWF valence and 
staggered (MILC) sea calculation. RBC-UKQCD's programme had previously published results 
from a single lattice spacing [6, II], and this result had been combined by Lellouch with the 
2-1- If staggered result of Gamiz et al. [51] in his world average. I choose not to include this 
staggered result in an average in light of the quenched discrepancy above. The RBC calculation 
had a preliminary update (since the Kaon conference) at Lattice 2009 with a second lattice spacing 
and a preliminary joint chiral and continuum extrapolation. These results are sufficiently significant 
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that they merit inclusion in these proceedings. This year Aubin, Laiho and Van de Water [52] used 
a mixed action approach with valence domain wall fermions on 2+lf of staggered sea quarks. 

A detailed comparison of both calculations is interesting in that estimates of systematic errors 
are in fact the dominant source of uncertainty. The statistical eiTors in both calculations are sublead- 
ing and it is worth emphasizing the subjective nature of the systematic enor estimates that dominate 
the quoted errors. The most important differences arise from the estimate of the renormaUsation 
and chiral-continuum extrapolation errors. 

ALV base their error on the non-perturbative renormalisation on the difference from mean 
field improved 1-loop lattice perturbation theory. However, RBC-UKQCD now use multiple non- 
exceptional momentum renormalisation schemes to gain control over their error estimate. 

Given figure 1, RBC-UKQCD estimated its chiral-continuum extrapolation error based on the 
difference between linear and SU(2) NLO extrapolation, figure 5. ALV apply SU(3), which is 
in principle less convergent than SU(2), and estimate the corresponding error based on varying 
a subset of analytic NNLO terms included in the fit. They do not include taste violation effects 
beyond NLO or non-analytic terms beyond NLO. It is likely that the difference here lies in the 
approaches to error estimation. 

The RBC-UKQCD result, however, is a preliminary conference submission and has not yet 
been given an estimate of finite volume effects, nor been published as even a preprint. Until these 
systematic errors have been published in peer reviewed journals 1 recommend continued use of [1 1] 

B^{n / = 2 + 1 , 2GeV) = 0.524(30) 





ALV 


RBC-UKQCD(Lat09) 


RBC-UKQCD(2007) 


6^S(«/ = 2 + l,2GeV) 


0.527(6)(20) 


0.537(19) 


0.524(10)(28) 


NP renormalisation 


3.3% 


2.4% 


2% 


Chiral/cont extrapolation 


1.9% 


2.4% 


2%, 4% 


statistical 


1.2% 


1.1% 


1.9% 


Scale & masses 


0.8% 


0% 


1% 


Finite volume 


0.6% 




1% 


total (quad) 


4.0% 


3.5% 


5.7% 



Table 6: Leading calculations of Bk with 2+1 f. Systematic error estimates are dominant, and a detailed 
breakdown of the sources error useful. RBC-UKQCD use a more convergent SU(2) approach to while 
ALV absorb taste mixings of heavy pseudoscalars with SU(3) ■ 

4.4 Renormalisation of Bk 

The dominant systematic error in Bk has become the QCD renormalisation. The four recent 
high quality 2f and 2-i-lf simulations [49, 50, 11, 52, 18] all make use of the Rome-Southamptom 
Rl-mom scheme to match between lattice and the MS regularisation schemes via an intermediate 
physically defined MOM scheme. 

By historical accident, the particular kinematic point selected was an exceptional momentum 
point, and suffered avoidably from [53] strong non-perturbative contamination. This paper demon- 
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strated that a better choice is to use non-exceptional momenta as the intermediate MOM scheme, 
and preliminary results have been presented at the lattice conference using a 1-loop calculation by 
Sachrajda and Sturm. In fact, multiple ostensibly equivalent 1-loop schemes were used and the 
spread gave a much more robust estimate of the systematic error at around 2.4%, leaving it as the 
(jointly) dominant source of error. 

While other components of the lattice determination of Bk will continue the rapid progress 
seen recently, the determination of this matching to two loops at this non-exceptional momentum 
point is required for reduction of the lattice error on Bk to below around 2%. 

The non-exceptional renormalisation point is also even more important to the robust determi- 
nation of the matrix elements of supersymmetry induced four quark operators [34]. 

5. Conclusions 

These are truly exciting times for lattice QCD as simulations push steadily closer to the physi- 
cal point in a large volume and in the continuum limit. There has been tremendous recent progress 
in a number of important areas; understanding of the chiral expansion of QCD has been greatly 
helped by Lattice QCD; Lattice QCD input to Vus is hugely important and the kaon bag parameter 
has been substantially better determined. 

The kaon is the ideal system for Lattice QCD; it contains rich and non-trivial CP violating 
physics at energy scales that can be simulated without large cut-off effects. Nascent work on elec- 
tromagnetic and iso-spin breaking effects was beyond the scope of this review; however, it appears 
that there are no barriers to the continued improvement in the precision of lattice simulations be- 
yond those posed by computing power. 
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